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ABSTRACT 

Bet«treen July 1973 and January 1974, three time 
transfer experiments using the Timation II satel- 
lite were conducted by the Division of National 
Mapping, A.C.T., Australia and the US Naval 
Research Laboratory, Washington, D.C., to measure 
time differences between the US Naval Observatory 
and Australia. Statistical tests showed that the 
results are vmaffected by the satellite's position 
with respect to the sunrise/sunset line or by its 
closest approach azimuth at the Australian station. 
Further tests revealed that forward predictions of 
time scale differences, based on the measurem.ents, 
can be made with high confidence. 

Measurements Against the Satellite Clock 

The results of the first two time transfer experiments 
between NRL and Australia have already been presented 
(Easton, Smith and Morgan, 1973, 1974) . 

The first series of statistical tests examined the residuals 
from a quadratic fit of the results TII-AUST, where Til re- 
presents the satellite on-board oscillator and AUST represents 
the local Australian time standard, in this case the National 
Mapping portable cesium standard DNM590 whose performance was 
linear with respect to a four component mean time scale. The 
measurements in July showed that the on-board oscillator had 
a constant aging rate during the experiment, and the residuals 
were normally distributed, while during the September run the 
residuals were not noraally distributed. The oscillator was 
evidently much less stable in the January 1974 run, and a 
simple curve could not be fitted. 

The curves fitted were: 

July rvin: TII-AUST = 5.371 + 4.7616{t-t) + 0. 25214 ( t-t) ^ 

microseconds, 

where t was the day of the year and t was the midpoint of the 
run (199.153). The standard deviation of the 38 residuals 
was 0.585 microseconds. 
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September run: TII-AUST =_38.940 + 3,9779(t-t) - 

0 . 00956 (t-t ) 2 microseconds, where t was 
265.544. The standard deviation of the 42 
residuals was 1.653 microseconds. 

For each run, the residuals were analyzed by a two-way 
analysis of variance with unequal cell sizes (Hamilton, 

1964) , the classifications being: 

(1) Effect of sunlight, to see if it affected either the 
satellite Clodk through temperature variation, or the signal 
travel time. 

This classification was stib-divided by the time of closest 
approach: 

(i) more than 2 hours before svmrise or more than 2 
hours after sunset; 

(ii) 1 to 2 hours before sunrise or 1 to 2 hours after 
sunset; 

(iii) O to 1 hours before stmrise or 0 to 1 hotirs after 
sunset ; 

(iv) o to 1 hours after sunrise or O to 1 hours before 
sunset; 

(v) 1 to 2 hours after sunrise or 1 to 2 hovirs before 

sunset; 

(vi) more than 2 hours after sunrise or more than 2 
hours before sunset. 

(2) Effect of azimuth, to see if the local surroundings 
had any effect. This classification was s\ib-divided by the 
azimuth at closest approach: 

(i) in the quadrant North to East; 

(ii) in the quadrant East to South; 

(iii) in the quadrant South to West; 

(iv) in the quadrant West to North. 

The number of observations falling into each cell, and the 
row and column means, are given in Tables I and II. The 
residuals from the fitted curves are shown in Figure 1. 

The analyses of variance presented show that there is no 
statistical evidence for effects due to the amount of sun 
illumination on the satellite-station path, or on the quad- 
rant of observation; nor is there any significant interaction 
between these two factors. No data were available to the 
authors for testing more precisely the effects of temperature 
and proton bombardment on the satellite oscillator. 
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January 1974 Results 

A third run was conducted in Januaary 1974, principally to 
provide an interpolating line for calibrating the Timation 
system against a USNO flying clock which visited National 
Mapping on 7 Decoonber 1973. On this run, oscillator in- 
stability, the restriction to one transmission frequency, 
and turning the oscillator on and off during the run de- 
graded the results of TII-AUST to the extent that compari- 
sons with the satellite clock could not be statistically 
analyzed. The time transfer comparisons with USNO were also 
degraded, as can be seen from Figure 2 , but the mean of the 
observations still proved useful. The results of fitting 
straight lines to the January measurements of USNO-AUST, 
together with combinations of all three rxins, are given in 
Table III. 

Inter-run Statistics 


Three series of tests were conducted to establish whether 
the runs, and their combinations, were statistically equiva- 
lent. 

2 

In the following descriptions, the c are variances of re- 
siduals after fitting straight lines of form 

USNO - AUST = a + P(t-t), 

2 

and the Cq are the variances of the observed rates P. The 
nxmiber of ^observations in each data set is denoted by n. 

(1) The first series compared the results obtained, on 
the one hand, by subtracting direct Australian observations 
from points interpolated between NRL observations, and on 
the other hand, by subtracting interpolated Australian 
observations from direct NRL observations. The tests were; 


(i) 


Equivalence of Sample Populations, i.e. whether 
NRL- interpolated samples were drawn from the same 
population as AUST- interpolated samples. 


Null hypothesis 
Alternative 


Test statistic 




f 


= a^/a 


2 

2 


(Fisher's F) 


Evaluation 


Accept if 0.60 < f < 1.67. 
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The f-column of Table IV shows that the populations were 
statistically equivalent at the 95% level for all except the 
September runs, which were nearly equivalent. 

(ii) Non-zero Significance of Rates, i.e. whether the 

rates of each run were statistically equal to zero. 

Null hypothesis : P = 0 

Alternative : P 7 ^ 0 

Test statistic ; t = CP]/(vp)^ (Student's t) 

: Accept if t < 1.96. 

The t-columns of Table IV show that the rates were very 
different from zero in all runs except January. This is due 
in part to the large standard error and small data set; but 
reference to Table III shows that the rate was indeed small, 
which is possibly explainable by the vagaries of the satel- 
lite oscillator which made the interpolation scheme unstable. 
It will be shown in a later test that the January rates were 
different from the rates determined from the other runs. 

(iii) Equivalence of Rates, i.e. whether rates obtained 
by interpolating NRL observations equalled AUST- 
interpolated rates. 

Null hypothesis Hq : P^^ ~ ^2 
Alternative ^2 

Test statistic : T - (P^^ - ^2^^^ (Student's t) 

where S = {[(nj^-2)CT^ + (n 2 - 2 )a 2 ^Xv|j^/cr^ + 

Evaluation : Accept H if T < 1.96. 

o 

The T-coliann of Table IV shows that the two interpolation 
schemes gave the same rates. 

(2) The second series of tests evaluated whether the 
midpoint of one run (a at time t in Table III) coincided 
with the value at t on the line fitted through another run, 
i.e. whether runs gave consistent values when extrapolated. 
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Null hypothesis: 
Alternative : 

Test statistic : 

d-2 degrees of 
freedom, where 


“l ” ^2 

®1 ^ °^2 ^ 2 ^^”^^ 
t = - [a 2 +p 2 with 

S = 


d = S^/la^/n^itij^+l) + CT2/n2(n2+l) 3 . 


This statistic is approximately distributed as Student's t - 
it is the incomplete Pisher-Behrens statistic (Welch, 1937; 
Hamilton 1964) . 


Evaluation : Accept if t < 1.96. 

Table V shows that forward extrapolation of the NRL- inter- 
polated samples is valid - even extrapolating from the July 
run into January is satisfactory at the 2 % level. The fact 
that the AUST-interpolated samples do not give such good 
extrapolation characteristics is attributed to the sparser 
Australian data sets - otherwise it is a little puzzling. 

(3) The third series tested the hypotheses that, for the 
NRL- interpolated samples, the sample populations and rates 
frcan each run and combination were equivalent; and similarly 
for the AUST-interpolated samples. All single runs satisfied 
a goodness-of-fit test for normality. 


(i) Equivalence of population variances. 


Null hypothesis 

Alternative 

Test statistic 
Evaluation 



^2 . „2 

^ ”2 

Accept if 0.60 < f < 1.67. 


The f-column of Table VI shows clearly that the July run had 
significantly Icwer variance than any other run or combina- 
tion, but that the other runs were, by and large, from the 
same population. 
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(ii) Equality of rates. 

Null hypothesis 

Alternative ^1 * ^1 “ ^2 

Test Statistic : T = [p^^ - (Student's t) 

where S = {[(nj^-2)a^ + in2~2)<3^io'^/(3^ + Vtnj^+n2-4]}^ 

1 2 

Evaluation : Accept if t < 1.96. 

The T-colxmin of Table VI shows that, for the NRL- interpolated 
samples, the rate determined from the January run was statis- 
tically different from the rates determined from all other 
runs and combinations, which were in turn statistically equal 
to each other. This confirms the result foiand in test (ii) 
of the first series. The poorer results obtained frcmi the 
AUST-interpolated samples confiann the second series tests 
wherein extrapolation between some run combinations was not 
valid. 

Comparison of Time Scales 

The tests described above all used a single cesium standard, 
DNM590, for the time scale denoted AUST. To demonstrate 
that the out-of-character results of the January anin were not 
due to a change of rate in this clock, a special Australian 
artificial time scale (AATS) was constructed, comprising the 
four cesitjm standards DNM590, the original Mount Stronlo 
standard DNM205, the newer standard NSL338 of the National 
Standards Laboratories, CSIRO, Sydney, and standard HP052 
maintained by Hewlett Packard (Australia) Limited, Melbourne. 
These clocks are all compared daily by ABC television com- 
parisons (Miller 1970) and were not stopped or adjusted in 
the period between 8 February 1973 and 23 May 1974. No other 
cesium standard in Australia satisfied both these conditions. 
The time scale was a simple unweighted mean of the four 
clocks, offset (in phase only) so that it agreed approximately 
with UTC (USNO) determined by flying clocks. 


An extrapolating ephemeris for AATS was constructed, using 
least squares straight line fits, in which: 


eCusno-aats] = E 


i=l 


and 


eC USNO-Clock^ 3/4 
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ECUSNO-Clock^] = e[uSNO-DNM 590] (by Timation) +ECDNM590-Clock^3 
(by Television) . 

Selected points on the graphs of the clocks against the arti- 
ficial time scale are shown in Figure 3. It can be seen that 
no significant rate change occurred in DNM590. Assuming that 
no such rate change occurred in OTC (USNO), the poor results 
in rate from the January run would reflect deficiencies in 
the Timation II technique* 

Table III includes the result of a visit by USNO flying clock 
PC572 in December 1973, and, in the colxamn headed USNO-AUST, 
gives the values obtained by inserting t = 7 December 1973 
in the various formulae obtained for different Timation II 
runs and combinations. For each interpolating system, the 
July-September combination gave the best agreement, which was 
0.31 microseconds for NRL-interpolations and 0.17 microseconds 
for HAST- interpolations. The 95% confidence interval at 7 
December for the NRL-interpolated July-September combination 
was +0.27 microseconds, and for the AUST-interpolated combina- 
tion +0.20 microseconds. Thus, on the ass;jmption that no 
error was attached to the flying clock result, the former 
set gave almost statistically' correct results, while the 
latter set showed excellent agreement. When the quoted error 
of +0,2 microseconds from the USNO certification was taken 
into accoxant, the NRL- interpolated result also became accept- 
able. 

Extrapolating the July-September combinations backwards to 
the date of the previous USNO flying clock visit on 8 February 
1973, the differences were 3.03 microseconds with +0.05 micro- 
seconds 95% confidence interval for USNO-interpolates, and 
3.12 microseconds with 0.37 microseconds 95% confidence in- 
terval for AUST-interpolates . When the errors of the flying 
clock measurements are taken into account, these results are 
not unsatisfactory. On the other hand, the agresinent between 
USNO-AATS by the extrapolating ephemeris and USNO-AATS by the 
flying clock measurement was 0.8 microseconds, which is regard- 
ed as very satisfactory over such an interval. It is unfortun- 
ate that no sxabsequent definitive flying clock trip has been 
made, as a further test of the predictive power of Timation 
would have been very beneficial, especially as time keeping 
in Australia during 1974 has been plagued with breakdowns 
both in a number of cesium standards and in the television 
network system. 

On the basis of the consistency of the NRL-interpolated 
July-September combination in both the extrapolative and 
interpolative senses, its agreement with USNO flying clock 
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measurements, and its consistency with a selected Australian 
time scale, the formula 4 in Table III was adopted as the 
definitive comparison between Australian clocks and UTC(USNO 
MC) , and was accordingly made the sole interpolating link 
between UTC(USNO MC) and the regular television-compared 
Australian mean time scale UTC(Aus). 

Conclusions 

The foregoing statistical analysis shows quite clearly the 
value of the Timation satellite for intercontinental time 
transfer at the silb-microsecond level. The major areas re- 
quiring particular attention are: 

(i) Long runs are required to establish rates reliably. 
The durations of the three runs were thirteen, 
seventeen and twelve days in July, September and 
January respectively, and even then the sparser 
results in January produced anomalies. 

(ii) The stability of the satellite oscillator must be 
good enough to carry interpolations over several 
hours. The poor January results establish this 
point forcibly. 

(iii) The superior results from the July run show that 
dual frequency transmissions (150 MHz and 400 MHz) 
do indeed reduce errors. 

(iv) Every effort should be made to have the data set 
as dense as possible. 

(v) There are factors affecting the stability of the 
oscillator which we do not yet understand, since 
an analysis of variance failed to reveal two 
possible causative, or perhaps correlated, effects. 
It is significant here that the residuals AUST-TII 
were not nomally distributed yet the residuals 
from USNO-AUST = (USNO-TII) - (AUST-TII) were norm- 
ally distributed, thus indicating the presence of 
a perturbing influence in the region of the satel- 
lite. No data was available to test accurately 
the hypothesis that the temperature aroxmd the 
crystal caused it to fluctuate. 

The Timation II results presented here have been incorporated 
into UTC ( AUS ) , so that predictions of the relationship be- 
tween Australian clocks and USNO can be made with confidence 
at the microsecond level. It is hoped that regular 
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observations of Timation III can be carried out - its im- 
proved clock should improve the statistics considerably and 
enable our geographically isolated clocks to contribute to 
International Atomic Time. 
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Figure 1: Hesiduals from Quadratic Pits of TII-AUST (NHL-Interpolated) , July and September 1973* 

The Horizontal Scale is Marked in UT Bs^s* The Vertical Scale is Marked in Mioroseconds 
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Figure 2i USNO-Aust (NRL-Interpolated). 

The Full Line is Line 6 of Table III, i*e* Combination of All Three Runs. 

The Bashed Lines are Lines 1,2, of Table III, i«e« Individual Huns. 

The Horizontal Scale is Marked in UT Bays. The Vertical Scale is Marked in Microseconds 
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Figure 3s Australian Artificial Mean Time Scale (AATS)-Clock for Each Component Clock. 
The Vertical Scale is Harked in Microseconds | its Origin is Arbitrary. 
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Analysis of Varianca» Aug/Sept Run 
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Validity of Extrapolating Between 



ifypothesis accepted at 256 level 
Hypothesis rejected. 













joValuation of Inter-ftun Consistency 
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Ifypothesis accepted at 2^ level 
Hypothesis rejected* 


